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Abstract 

We explore the nonlinear classical dynamics of the three-dimensional the- 
ory of "New Massive Gravity" proposed by Bergshoeff, Hohm and Townsend. 
We find that the theory passes remarkably highly nontrivial consistency checks 
at the nonlinear level. In particular, we show that: (1) In the decoupling limit 
of the theory, the interactions of the helicity-0 mode are described by a single 
cubic term - the so-called cubic Galileon - previously found in the context of 
the DGP model and in certain 4D massive gravities. (2) The conformal mode 
of the metric coincides with the helicity-0 mode in the decoupling limit. Away 
from this limit the nonlinear dynamics of the former is described by a certain 
generalization of Galileon interactions, which like the Galileons themselves 
have a well-posed Cauchy problem. (3) We give a non-perturbative argument 
based on the presence of additional symmetries that the full theory does not 
lead to any extra degrees of freedom, suggesting that a 3D analog of the 4D 
Boulware-Deser ghost is not present in this theory. Last but not least, we 
generalize "New Massive Gravity" and construct a class of 3D cubic order 
massive models that retain the above properties. 



1 Introduction and Summary 



In this work we study massive gravity in 2+1 spacetime dimensions. A well-known 
theory of this is Topological Massive Gravity [1]. Here, we focus on the nonlinear 
dynamics of "New Massive Gravity" (NMG) proposed by Bergshoeff, Hohm, and 
Townsend, [2]. Similar studies on the nonlinear structure of massive gravity in 3+1 
dimensions have by now reached a certain level of maturity. Since these develop- 
ments are relevant for the motivation and content of the present paper, we begin 
with a brief outline of the most important insights gained. 

Fierz and Pauli (FP) [3] constructed a ghost-free and tachyon-less linear theory 
of massive spin-2 in Minkowski space, describing 5 degrees of freedom in 4D (2 
degrees of freedom in 3D), that is consistent with the corresponding representation 
of the Poincare group. The FP theory has no continuous limit to the massless theory 
[1], while the continuity can be restored in its nonlinear extensions [5]. However, 
a generic nonlinear extension in 4D suffers from the so-called Boulware - Deser 
(BD) ghost: the Hamiltonian constraint, that would restrict the number of degrees 
of freedom to no more than 5 for any background, is lost at the nonlinear level 
[6]. The consequences of this can be transparently seen in the decoupling limit 
[7]; here, the problem emerges in the equation of motion for the helicity-0 mode 
that exhibits ill-posed Cauchy problem due to nonlinear terms with more than two 
time derivatives acting on a single field [HI |9]. A toy example of this in 4D is the 
Lagrangian 

where, A5 is a certain scale (composed of the graviton mass and Planck mass), 
and T denotes the trace of the external energy-momentum tensor. The equations 
of motion that follow from ([1]) involve nonlinear terms with more than two time 
derivatives acting on vr. This leads to an additional (sixth) mode with negative 
kinetic term propagating on nontrivial backgrounds for the helicity-0 mode (e.g., 
that of a collapsing spherical source) [TUl El IS]; this mode is identified with the BD 
ghost. This follows directly from the standard Ostrogradski instability arguments 
and can be made manifest here by noting that the above theory can be equivalently 
represented by the following lagrangian , 

C = i« - - eJ^l^'^Ai/^ + - -1-,PT. (2) 

The latter lagrangian is obtained from ([1]) in two steps. First one introduces an 
auxiliary field A and a corresponding lagrangian C, which involves no more than 
single derivatives per field. The Lagrangian C is defined so that integrating out A, 
one recovers the original vr-lagrangian ([T]). On the other hand, making the further 
field redefinitions \^ = ip and n = (p — ip brings C to the form, given in Here 
e = ±1, depending on the sign of □(</) — ip). 
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It is manifest from (|2]) that the additional degree of freedom ip is inevitably a 
ghost; however, its mass is set by {ip)~^^'^, and although this mass is infinite when 
expanding around a Minkowski background, which explains why its presence was 
not seen in the original FP analysis, it can drop down to physically accessible scales 
when expanding around a non-trivial background, [8|. 

The simplest example of a nonlinear theory for a helicity-0 mode graviton that 
does not have the above problem was found in the context of the DGP model pjj 
in Ref. pS]- Its Lagrangian reads as follows: 

1 (dnYDn 1 

The specific structure of the cubic term prevents the appearance of more than two 
time derivatives in the equations of motion for vr, thereby rendering it free of the sec- 
ond, ghost-like, degree of freedom on any weak asymptotically flat background, [12] . 
The Lagrangian ([3]) has been generalized to incorporate the quartic and quintic self 
interactions with similar properties - the so-called Galileon terms, [13] . These terms 
have been shown to naturally arise in probe brane setups [Hj and their extensions 
to multi- Galileon theories have been studied in [151 [IHl ITT] . 

More importantly for the present work, Refs. [T8l [19] showed that a certain class 
of theories of 4D massive gravity, contain all the Galileon terms in the decoupling 
limit. In these theories, an infinite number of nonlinear terms can be resummed 
and a Lagrangian containing just a few terms can be obtained [20j. Hence, the 
BD problem can be addressed in the full-fledged Hamiltonian approach. So far the 
Hamiltonian has been calculated up to and including quartic order in nonlinearities 
and beyond the decoupling limit, showing that the BD problem does not arise at 
this level ^0. 

The purpose of the present work is to both study the nonlinear interactions of 
the helicity-0 mode in the decoupling limit of NMG [2], determining whether or not 
it gives rise to the problematic nonlinear terms, as well as to address the issue of 
the BD ghost in the full theory, away from the decoupling limit. Since NMG is 
a 3D theory, one would expect to be able to do more than in 4D, as well as gain 
some additional intuition about the 4D constructions. The divergences arising in 
3D being less severe, NMG can give some valuable insight on the renormalization 
of gravity, [22]. Furthermore, a DBI extension was proposed in Ref. [23], and a 
remarkable connection with AdS/CFT was established, [23]. The existence of AdSs 
Black Hole solutions also makes this theory especially interesting, |25j . 

In this work, we will focus instead on the stability of the theory and show that 
the interaction Lagrangian for the helicity-0 in the decoupling limit of NMG reduces 

^ A potentially dangerous term found in the quartic order of the Hamiltonian construction in 
Ref. [5] has a very special structure [2D], and due to this, can be removed to higher orders by a 
consistent nonlinear field redefinition, as shown explicitly in [20j. The terms that could not have 
been removed, do cancel out automatically in the quartic order |20] . We also note that the term 
of concern of Ref. [21) . found in that work in the Lagrangian formalism, is nothing but the term 
already found in f8^ in the Hamiltonian formalism, that was addressed above. 
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to the three-dimensional version of Lagrangian ([3]), supplemented by an equation 
which determines the tensor modes. Hence, the theory has a well posed Cauchy 
problem. The nonlinear term in seems to arise universally: It appears in the 
context of the 5D DGP model, as found in [12], in 4D massive gravity [T8l IT9| 120] , 
and now in 3D NMG. In what follows we demonstrate that in NMG, the helicity-0 
mode coincides with the conformal mode of the metric in the decoupling limit. This 
motivates us to study the dynamics of the conformal mode in the full theory. In 
ordinary OR, it is well known that the conformal mode has the wrong-sign kinetic 
term. This is not a problem in GR since this mode is not dynamical (it can be 
removed via the gauge freedom and constraints). In NMG the conformal mode is 
dynamical, but the GR term itself has a wrong sigiJl, hence, the conformal mode 
obtains a correct sign kinetic term. Furthermore, the nonlinear terms of NMG give 
rise to interactions for the conformal mode that are straightforward generalizations 
of the Galileon interactions, implying that they have a well-defined Cauchy problem. 
Moreover, we give a non-perturbative counting of degrees of freedom that demon- 
strates that in NMG, even away from the decoupling limit or from conformally flat 
metric configurations, there are no extra modes beyond the two of a massive 3D 
graviton. This confirms that the 3D analogue of the BD ghost does not arise nonlin- 
early. Last but not least, we present a class of 3D cubic order theories that generalize 
the NMG construction. 

The paper is organized as follows. Section 2 serves as a review where we present 
some well-known properties of massive gravity in 4D, adapted to three spacetime 
dimensions. We illustrate that all known potential problems that arise when at- 
tempting to give a 4D graviton mass, persist in generic three dimensional models as 
well. The reader who is well acquainted with the technology and issues involved is 
encouraged to skip directly to section 3 where we investigate the nonlinear dynamics 
of NMG. We start by studying the decoupling limit of the theory, and arrive at a 
well-defined, ghost-free theory. We then move to the dynamics of the conformal 
mode, and find that it coincides with the helicity-0 graviton in the decoupling limit. 
We show next that for conformally fiat metrics the full theory is free of any ghosts. 
Moreover, we give the non-perturbative arguments in favor of no-ghost propagation 
beyond any limit or approximation. In section 4 we look at NMG from the per- 
spective of a generalization of the linear FP model, demonstrating how exactly the 
cancellation of the BD ghost happens at the cubic level. Finally, we construct a 
class of generalizations of the cubic theory, which exhibit similar properties. 

2 Ghosts and Strong Coupling in Massive Gravity 

In this section we review some of the known results from massive gravity in 4D and 
introduce the formalism used in this paper. We begin by analyzing the Fierz-Pauli 
model at the linear level. We then extend the discussion to include nonlinear terms. 

•^This is acceptable as there are no propagating degrees of freedom in 3D massless GR. 
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2.1 Linear analysis 

The FP model is the unique theory at the hnearized level, free of ghosts or tachyons 
[26] and propagating only the two degrees of freedom of a massive spin-2 particle 
in 3D. In terms of the metric perturbation h^y = g^i, — rj^,^, the FP Lagrangian is 
given by, 

C = Auph^^Sphp^ - ^Mpm\h''''h,, - h') + h'^'^T^,, (4) 

where m sets the graviton mass, T^i, denotes an external energy-momentum source, 
and Sjl'^ is the linear Einstein operator 

^^"V = -^(DV - d.^d'^h^, - d.d'^K^ + d,Ah - r],,,Uh + r^^.d^dph''^), (5) 

with all indices contracted with the flat metric. 

In order to single out the propagating helicity-1 and helicity-0 modes of the 
massive graviton, the metric perturbation can be conveniently decomposed at the 
linear order in the following way, 

where h^,^, and vr encode the (canonically normalized) helicity-2, 1 and compo- 
nents, respectively (note that in 3D the helicity-2 mode does not propagate). With 
this field content, the theory is invariant under the usual linearized diffeomorphisms 

h^,u h^y + df,^u + dui^, Af,-m^^ , (7) 

as well as an additional abelian U(l) symmetry, under which the vector and scalar 
modes transform as, 

A^,-^ Af, + df,C, TT ^ TT - m C- (8) 

The latter invariance, which at this stage is introduced somewhat artificially, will 
turn out to provide a very convenient bookkeeping tool when studying different 
decoupling limits of nonlinear massive gravity. 

As remarked in the introduction, the linear FP theory does not possess a con- 
tinuous massless limit. To see this, we note that in terms of the different helicity 
components, the m — )■ limit of the theory is given (up to a total derivative) by the 
following expression, 

Cm=o = —h^^'S'.lhp. - If^'F,. - h'^'^id.d.Tr - r/^^nvr) + -^h^^T^,, (9) 

where F^^, denotes the usual abelian field strength for A^. The mixing between the 
tensor and scalar modes in ([9]) can be eliminated by a linear conformal redefinition 

= hf,u + 2r]^^7r, (10) 
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which brings (Q to the following form 

c^o = -ift'-f-zv, - If'-f,. + 2.n. + ^w-r,,. + ^.r. (ii) 

It is clear from the latter Lagrangian that the helicity-0 part of the massive graviton 
does not decouple from matter even in the m — )■ limit, leading to the famous vDVZ 
discontinuity [?]. This is an 0{1) modification of the gravitational interactions as 
compared to GR, at least in the regime of validity of the linear approximation. 
However, as first pointed out in [5], in nonlinear theories of massive gravity this 
approximation typically breaks down at a parametrically large distance ry from 
localized sources, allowing for their phenomenological viability. In fact, it is the 
nonlinear dynamics of the scalar mode itself that screens its contribution to the 
gravitational potential within the Vainshtein radius, ry, restoring agreement with 
GR. In a generic nonlinear extension of the FP theory, however, the same nonlinear 
self-interactions of vr are responsible for a number of theoretical problems, such 
as the propagation of ghosts and the ill-posedness of the Cauchy problem, strong 
coupling at parametrically low scales, and etc. Before turning to the discussion of 
theories that avoid these problems, it is instructive to have a closer look at their 
origin by considering the minimal nonlinear extension of the FP model, to which we 
turn next. 



2.2 Fierz-Pauli model at the nonlinear level 

Let us "naively" continue the theory of the previous subsection at the nonlinear 
level by keeping the graviton potential unchanged, while completing the derivative 
self-interactions to the full Einstein-Hilbert action, 

C = M^^gR - ^ [hl^ - h^) . (12) 

Here index contractions on the metric perturbation are performed using the flat 
metric rj^'^ . In order to single out the high-energy degrees of freedom in f|T2l) . we 
employ a trick a la Stiickelberg [27]: For the purpose of dealing with the theory at 
cubic order, it is more convenient to restore general covariance by introducing the 
decomposition for h^^, patterned after the nonlinear gauge transformation of the 
metric perturbation 

VWm^ ^ Mpm' ' 
where h^y and are further decomposed as 

d V^h + d V^h + V^d h 
V = V+ " TMpt^ ' V, = mA, + c?,vr. (14) 
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The fields h^^ = h^i, — 27] ^^n, A^^ and vr will describe the canonically normalized 
tensor, vector, and scalar modes at high energies, respectively. Plugging the decom- 
position ( I14p into the Lagrangian ( I12p . and considering the limit 

Mp ^ oo, m ^ 0, (Mp/^m^)2/9 = A9/2 = finite, (15) 

one recovers a particular high-energy ( "decoupling" ) limit of the theory. Performing 
the conformal shift ( TTOl) of the helicity-2 mode, (fT2|) reduces up to a total derivative 
to the following expression, 

Cdec = -h'-'S'.fh,^ - If^'F,, + 27rn7r + {{Dnf - (□7r)(9^9,7r)2) 

+ I Ji^^T.^ + -^=7iT. (16) 

As noted in the introduction, the cubic decoupling limit of the FP theory is amended 
by the helicity-0 self-interactions of the form {d'^ir)^. They successfully screen the 
scalar contribution to the gravitational potential of a point source of mass M inside 
the Vainshtein radius ry , given by 

which restores the agreement with General Relativity. This however happens at 
the expense of introducing a ghost in the theory: Although infinitely heavy on a 
fiat background, the mass of the ghost becomes low enough to be disruptive around 
any reasonably localized source [TJ [H |9]I^ . Even if one does not consider an external 
source, the cubic self-interactions of the scalar field, containing more than two time 
derivatives, cause the Cauchy problem to be ill-posed for the theory at hand. 

Recently, a class of theories of (four-dimensional) massive gravity that avoid these 
problems have been constructed fl8\ [T9| [20] . These theories modify the graviton 
potential order-by order, so as to cancel all potentially dangerous self-interactions of 
the scalar field. As a result, one obtains a sensible effective field theory for a massive 
graviton, with the cutoff given by A3 = {Mprn^Y^^ in four spacetime dimensions. 
Moreover, the decoupling limit of the theory, obtained by keeping A3 finite while 
setting the graviton mass and Planck constant to zero and infinity respectively, is 
free of ghosts below the cutoff. Furthermore, the form of the decoupling limit is 
unique at the quartic order; in other words, any nonlinearities in the potential of 
order higher than the quartic one have no effect on ilQ. 



^Although the analysis in these references is performed in four spacetime dimensions, we con- 
firmed that the conclusions remain intact in (2-|-l)D as well. 

^Interestingly enough, a recently proposed nonlinear completion of Fierz-Pauli massive gravity 
[25] automatically produces exactly the right structure of the potential at the cubic level, so 
as to fall into the category of such models, |18j . 
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3 New Massive Gravity 



In this section we analyze the nonhnear dynamics of NMG. As we show at the end of 
the section, viewed as a generahzation of the FP model, NMG represents a curious 
example. In addition to introducing a potential for the graviton, it also modifies 
its derivative self-interactions beyond the ones originating from the Einstein-Hilbert 
action, leading to a well-defined decoupling limit similar to the four-dimensional 
case of [m [ini |20]. Most importantly, however, the full theory can be shown to be 
free of the BD ghost. We turn to the discussion of this remarkable property next. 
One way to define NMG is through the following action 



Snmg = Mp j d^x y/g 



:i8) 



Here G^y and R are the usual Einstein tensor and Ricci scalar for the metric g^j^y, 
respectively. All the indices are raised with the inverse metric g^^ and /^,^ represents 
an auxiliary covariant tensor field (/ = g^^f^u) that can be integrated out to yield 
a particular higher-derivative extension of the 'wrong sign' GR The equations 
of motion, obtained by varying ( TT51) with respect to g^i, and /^,^ (we keep the overall 
factors of which will be essential for the perturbative analysis of the equations 
beyond the linear level) are given by the following expressions 

VaG^'" = \vgg"'r^G^p - ^fl^.G'^'' - ^ZG^^' + ^v^(/i?^^ - rR) 

+\^^./gg"'{r^U - f) - l^'^if^r - F'f) , (19) 



VgG'"' + ^Vg^'ir - g'^f) = o, (20) 

and admit a Minkowski vacuum with gf^i, = rj^^ and f^y = (the Einstein and 
Ricci tensors are taken as functions of the metric gf^^ in the latter equations). From 
Eq. ( !T9|) we see that at the linear level the perturbations h^i, and f^u of the metric 
and the auxiliary field satisfy 

g2(/^-/) = o, 

where G^j^] denotes the linearized Einstein tensor. Since it does not propagate any 
degrees of freedom in three dimensions, the only solution to the latter equation 
is (up to a gauge, which we choose to be zero) h^^ = f^^. This, in combination 
with Eq. fl20|) . yields the usual FP equations of motion for f^^ in the linearized 
approximation. 

^Note the minus sign in front of the Einstein-Hilbert term in (1181) . which is characteristic of 
Topologically Massive Gravity as well 



7 



3.1 Exact decoupling limit of NMG 

Before turning to the full theory, we start by considering the decoupling limit of 
NMG, defined as, 

Mp cx), A5/2 = (A/iWpm^)^/^ = fixed. 

We consider the (a posteriori justified) ansatz 

"- = ;te' /.^ = ^ + v.vi + v„v;. (21) 

where is a vector field which will encode the helicity-0 and -1 degrees of freedom 
of the massive three-dimensional graviton, and V denotes the covariant derivative 
associated with the metric g^y. Plugging the latter decomposition into ( ITSl) . we 
obtain, up to a total derivative, the following expression for the Lagrangian of NMG, 

1 — 7?7^ — — — 



C-NMG — Mpy/g 



2 

m 



2 2 
TD — Tf) 



2 



(22) 



In here, the last term comes from the commutator of covariant derivatives and the 
quantities R and G^^^, as well as the covariant derivative itself, are associated with 
Qf^i^. Decomposing further into the canonically normalized helicity-1 and -0 modes 



and vr. 



V = \ f23) 



one arrives at the Lagrangian, which up to a total derivative can be written as, 
Cnmg = 



-MpR - ^/WprO^, - ^m\rf,, - f ) _ mr{V,A, 



+ W^A^A, + -R^^-V^nA, + ^R^^V^TiVy-K 
2 m 

It is straightforward to check that at the linearized level the theory propagates the 
two degrees of freedom (encoded in A^ and vr) of the massive 3D graviton. At this 
point we take the exact decoupling limit Mp — > 00 defined by keeping the scale A5/2 
fixed, giving rise to the following Lagrangian, 

+ 47^i?(i)^-9,7r9,7r, (24) 
^5/2 



where R^f!^} denotes the hnearized Ricci tensor for /i^i,, F^,^ is the usual abehan field 
strength for and all indices are contracted with the fiat metric. The auxiliary 
field fij,u, being a Lagrange multiplier in the limit at hand, imposes the constraint 

(Sh)^^ = -df.d^Tv + r^^^DTT, (25) 

which is solved by, 

h^u = 27r 1]^,^ + gauge transformation. (26) 

Plugging the latter expression back into the Lagrangian and extracting a total deriva- 
tive, one finally obtains the exact decoupling limit of the theory, 

C%la = 2'nU'n - \{d^YU'n - \fI^. (27) 

The decoupling limit of NMG therefore contains the helicity-0 mode with the cubic 
Galileon self-interaction, together with a free helicity-1 mode! The tensor modes 
hfj_y and f^^ have both disappeared in this hmit because they become massless, and 
by standard arguments massless spin 2 fields in 3D carry no propagating degrees 
of freedom (although in the presence of sources they will still contribute to the 
interaction energy). Since h^^ = 27ir]^^+ gauge transformation, we automatically 
generate a coupling to matter that is of the form oc vrT if external sources are 
considered. Remarkably, as we show below, this is an example of a theory of a 
massive 3D graviton with modified derivative self-interactions, free of ghosts to all 
orders in the decoupling limit. Even more importantly, we show next that the 
absence of the BD ghost persists even away from this limit. 

3.2 Dynamics of the conformal mode 

Before presenting arguments in support of this last statement, as a preliminary 
step it is very instructive to have a look at the dynamics of the conformal mode 
in NMG. Integrating out the auxiliary field f^y in f fTSj) . one arrives at the original 
representation of the NMG action |2], 



Cnmg — Mpy/g 



If 3 . 

R H t: R^'^Rau — -R 



(28) 



„ 1 -II, -11,^1/ 

Considering conformally fiat configurations of the metric 

9>.u = ^ = e-^/^, (29) 

the Lagrangian (1281) can be rewritten in terms of the field 0, as 
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where all indices are contracted with the flat metric. Curiously enough, the A5/2 de- 
coupling limit of this theory coincides in form with the decoupling limit Lagrangian 
of the helicity-0 mode in Eq. ( 12 7p 

£^^,G^20n0-i(90)^n0, (31) 

identifying the helicity-0 part of the massive graviton with the conformal mode in 
this limit. One can show that even away from the decoupling limit, the theory given 
by (|5UI) does not propagate ghosts as one would naively assume from the presence 
of higher derivative interactions. 

Indeed, up to a total derivative, Eq. (|30|) can be rewritten as. 

Naively, the second and third terms in the above Lagrangian, involving more than 
two time derivatives, might lead to the presence of ghosts (or equivalently, the ill- 
posedness of the Cauchy problem) in the theory. One can however show that the 
specific structure of these operators makes them harmless. The second term is 
conveniently expressed in terms of the Levi-Civita symbol as follows, 

oc e-^/^£^"'"^£,;3p<x d.d-'ct) d^d^<p. (33) 

The antisymmetry of the Levi-Civita symbol can then be used to show that no 
terms with more than two time derivatives are present in the equations of motion. 
Similarly, noticing that the third term in Eq. f l32p includes a factor of a peculiar 
form 

it becomes fairly straightforward to show that it does not cause the Cauchy problem 
to be ill-posed either. Indeed, this latter expression, being (up to a total derivative) 
equivalent to the DGP galileon □0(5^(/))^, is well known to lead to no more than two 
time derivatives per field in the equation of motion. As can be straightforwardly 
checked, the factor of e~'i'l^f^ in front does not alter this property. 

3.3 No ghosts in New Massive Gravity 

Finally, we give an argument demonstrating the absence of the BD ghost in the 
full NMG. The argument we shall give is based on counting non-perturbatively the 
degrees of freedom utilizing the symmetries. The key observation is that there is a 
non-perturbative generalization of the f/(l) symmetry described in Eq. ([8]) which 
combined with the existing 3D reparameterization invariance is sufficient to demon- 
strate that only the 2 propagating degrees of freedom remain non-perturbatively. 
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The Lagrangian fl22l) is invariant under the usual reparametrizations of coordi- 
nates (3D diffeomorphisms) , under which the metric and the auxihary field f^i, 
transform as tensors, while is a vector. One can use this symmetry to eliminate 
all potentially propagating degrees of freedom in the 3D metric Qf^u- By introducing 
the vector field we also introduce a new three-parameter local invariance (which 
we refer to as secondary linearized 3D diffeomorphisms) , under which the metric g^u 
remains unchanged, while f^i, and transform as, 

-> Uu + ^/Wpiy.i, + , v^-^v^-i,. (34) 

It will prove helpful to slightly rewrite the Lagrangian f l22|) in the following form. 



^NMG — \f9 

1 



-MpR - ^p~rG,, - \m\~r~f^, - p) - \f>^'F^. 



1 TD _ _ _ _ _ 1 Td'^ — — 

where the following set of field redefinitions has been used. 



(35) 



= ^/Wpmv^, f,, = f,, - ^-^(v;k - g.uV'^V^), (36) 

and F^y denotes the abelian field strength for V^. The symmetry (IMI) induces the 
corresponding transformation on f^^, that leaves the action invariant. One can use 
this freedom to ehminate all potentially propagating d.o.f's in f^y, leaving the 
only propagating field in the theory. 

Generically, would propagate three degrees of freedom in three dimensions; 
however, as we shall argue below, NMG is special in this sense, propagating only 
two degrees of freedom in the latter field. The only place in the Lagrangian f l35p 
where the time derivative of Vq appears, is the last term of the second line and the 
first term on the last line, which however include Vq only linearly. This means that 
the canonical momentum, conjugate to Vq is independent of Vq itself, enforcing the 
primary constraint on the system. This reduces the number of propagating degrees 
of freedom in to two, which, due to the arguments presented above, means that 
these are the only two dynamical degrees of freedom in the full theory of NMG. To 
make this argument clearer we can define 

V^ = A^ + d^n, (37) 

where is a vector under usual 3D diffeomorphisms, and transforms as Afj_ — )■ 
Afj, — C,fj. under the secondary linearized 3D diffeomorphisms, and tt is similarly a 
scalar under the former and invariant under the later. Since the Lagrangian (!35|l 
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depends only directly on V^, it is manifestly invariant under a new additional U{1) 
symmetry 

A^^A^ + d^C, vr ^ vr - C, (38) 

generalizing the result in the linearized theory Eq. ([H]). Since the above Lagrangian is 
only linear in Vq it is in turn guaranteed to be only linear in vf. Any such interaction 
will lead to second order equations of motion for vr. However, one may worry about 
defining a conjugate momentum for vr. That this can be done consistently, can be 
seen by noting that it is always possible to integrate by parts to put the action in 
a form where there are no more than single time derivatives acting on the fields 
in the Lagrangian. To see this explicitly we note that the only term in the above 
Lagrangian where the vf terms arise, is in 'gauged' Galileon interaction on the last 
line 

^0 = -\V9^^V''{V,V, + g.^V'^V^) , (39) 

as well as in 

^1 = -^/gr (v^K - ^7^.v"K) = {v,v,) ir - hn ■ m 

On integration by parts we see that the first term is equivalent to 

^0 = -lv9iV"V^)'^% = -\iV'^Va)d,i^Vn . (41) 

It is now easy to show that the would be problematic vr coming from the Vq part is 
a total derivative and consequently can be removed by an integration by parts 



1 



+ terms with no Vq, iV or N\ (42) 

where A^ is the ordinary lapse (7"° = —l/N"^, g^i = Ni the associated shift and 
9ij = lij- This is precisely for the same reason that the covariant DGP Galileon 
interaction □7r(97r)^ term gives well defined equations of motion. Indeed setting 
= we see easily that this term is equivalent to the DGP interaction for tt: 

(\/"\4)V'^\/^U,=o = nTx{dTxf . (43) 

Implicitly these integrations by parts determine the analogue of the Gibbons- Hawking 
boundary terms for NMG. A similar argument can be employed for £1, where all 
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potentially problematic terms can be seen to be a total derivative: 

Ci = -^N„ - V„JV - ViN' + Vi + . . .j + ■■■ 

= -V77«/«(5o(|)-Va(^))+-- (44) 

SO that this term is also free of any time derivative on the Lagrange multipliers, 
>Ci = -^0 {Vo - V,N')^ + terms with no N, N\ }oo or (45) 

Again the key feature to stress is that at the same time as removing the Vq terms 
we can also remove the N, TV*, /oo and foi terms ensuring that A^, A^', /oo and 
foi act as Lagrange multipliers for the normal 3D diffeomorphism Hamiltonian and 
momentum constraints. Similarly since there are by construction no Aq terms, then 
Aq acts as the Lagrange multiplier for the U{1) 'Gauss law' constraint. 

Putting this together, since all the double time derivative terms for tt can be 
removed by integration by parts, and by the same means the Lagrangian can be 

made independent of A", A^*, /oo, foi, ^o, we are guaranteed that the equations of 
motion are well defined, the 7 constraints arc preserved, and that there is a well- 
defined Hamiltonian formalism. The key point then is that following these steps the 
configuration space action can be written in a form in which it is a function of the 
6+6+3+1=16 variables g^u, f^u, A^,n each of which has well defined equations of 
motion. In addition there are 3+3+1=7 first class symmetries corresponding to 3D 
diffeomorphisms, secondary linearized 3D diffeomorphisms and the additional U{1). 
Associated with each symmetry there is a gauge freedom and a constraint which 
allows us to remove 2 degrees of freedom per symmetry. Thus the total number of 
degrees of freedom non-perturbatively is 16 - 7 -7 =2. Since 2 is the correct number 
of physical polarizations of a massive spin 2 field in 3D this confirms the absence of 
the BD ghost to all orders, consistent with our previous arguments which were only 
valid in certain regimes. 

4 Explicit Computation of the Decoupling Limit 
at Nonlinear Order 

In this section we analyze the decoupling limit of NMG explicitly to cubic order. 
We show that one indeed recovers a ghost-free theory, governed by the scale A5/2. 
We do so by exphcitly demonstrating how the terms corresponding to the A9/2 and 
A7/2 all vanish. We then generalize the results and construct a new class of theories 
which are ghost-free to cubic order. 
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4.1 NMG at Nonlinear Order 



To obtain the NMG Lagrangian to cubic order one must integrate out the met- 
ric perturbation h^y from the equations of motion of NMG f|T9l) - fl20l) beyond the 
quadratic order in the Lagrangian. This is a rather lengthy computation and it is 
outhned in the appendix. The resulting cubic effective Lagrangian for f^y can be 
written in the following form, 



4 



+ 7;flf-Vl-7;f)- (46) 



Here 7^,^ = rj^y + f^y is the analog of the full metric with f^y viewed as metric 
perturbation, and Mp{^R{'^)Y'^'> denotes the cubic in f^^ part of the corresponding 
Einstein-Hilbert action. At the quadratic level, we of course recover the FP action; 
a remarkable peculiarity of the cubic part, however, is that it modifies GR not just 
by the graviton's potential, but also by the deformation of the nonlinear Einstein- 
Hilbert part itself. 

According to the discussion of section 2, in a generic nonlinear extension of 
the FP action we should expect A9/2 to be the scale, governing the dynamics of 
the helicity-0 mode. On the other hand, we have shown in the previous section 
that A5/2 represents the lowest scale, by which the Galileon self- interactions of the 
scalar graviton are suppressed in NMG. The theory must therefore possess a special 
structure, which eliminates the scale A9/2, as well as the intermediate scale A7/2 = 
{yj Mpm^Y'f'^ ^ from its dynamics. To show this, let us plug the decomposition ( ITSll 
(with the obvious change h f, while keeping the same notation for the vector 
and scalar modes) into the cubic Lagrangian Eq. fH6|) . In the limit where all scales 
greater than A7/2 are sent to infinity, one recovers the following expression for the 
action, 

^9/2 

+ ( {d^'A''d''d,7i + d^APd^dpTi) + ^d^dPr^d^dpA S^y^Up. (47) 

\^7/2 ^9/2 / 

Here f = f — ^tci] represents the familiar conformally shifted tensor mode, so that 
in terms of the new fields the U{1) gauge symmetry of ([H]) is given by, 

Uu + ^m^Vt^y, A^^ Af, + df,^, TT ^TT-m^. (48) 
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Under these transformations, the Lagrangian (147|) is invariant up to terms that 
vanish in the decouphng hmit, 



The TT self-interactions in the second hne of Eq. fH7|) combine to a total deriva- 
tive and therefore drop out. Furthermore, extracting a total derivativ^ from the 
dAdnddTT-ty^ie operators and using the gauge invariance Eq. (HHl) to impose the 
condition d^A^ = 0, (HTl) can be rewritten as. 



h/2 



^/ = --rS''^i!U--F^''F,, + 2TiUTi + ^d,A,d^^d^TiUT, 



+ ( ^ {d^Afd^'dpTi + d^'Afd'-dp'K) + ^d^^dPixd^dpTA s^l^Up. (so) 

. ^7/2 ^9/2 / 

The remaining interactions can be removed by a field redefinition of the helicity-0 
and helicity-2 fields, 

7r^7r-A^d^A,d''d'7r;f^,-^f^, + -Ij^ {d^A^d^dpTT + d,APd^d,7r) 

4^7/2 ^7/2 



1 

^9/2 



+ -jgj^d^dPnd^dpTT, (51) 



As promised, all the dangerous interactions associated with the A9/2 and A7/2 
scales have dropped-out in this limit. This leaves A5/2 to control the nonlinear 
dynamics of NMG, as shown in subsection 3.1. 

It is important to note the significance of the U{1) gauge invariance in the pre- 
ceding discussion. In fact, precisely due to this freedom, any operator of the type 
dAd'^nd'^Ti at the scale A7/2 can be removed (up to a total derivative) by a redefi- 
nition of the helicity-0 graviton. This means that once the problematic (9^7r)^-type 
operators at the scale A9/2 are removed by a judicious choice of the coefficients of the 
different terms in the Lagrangian, the A7/2 scale becomes automatically removable, 
leading to the decoupling limit, governed by As/J^. 



^The corresponding total derivative is given by 
7 {df'A''d^d''Trd,daTT - d^A'^d^d^TTOTr ~ ^d^^A^id^d^Tr)^ + i9''A^(n^)2) /A^/^. 

^Mathematically, this traces back to the fact that the U{1) transformation of the scalar in ([8]) 
includes a factor of the mass m, connecting the scale Ag/2 to A7/2 = m/Ag/2. Because of this, 
tuning the theory so as to remove the former scale automatically results in the removal of the latter 
one, as was for instance encountered in 120) . 
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4.2 General Cubic-Order Ghost-Free Massive Gravity 

Finally, we extend NMG at cubic order considering a more general set of deriva- 
tive self-interactions of the graviton. The generic type of cubic Lagrangians we 
wish to consider is written in terms of the Einstein-Hilbert term (with possibly a 
non-standard coefficient) together with reparametrisation invariance-violating con- 
tributions involving the metric perturbation h^^i, = g^y — rj^^, 

,2 



C = Mp 



p{^Rp - ^ {hi, -h' + chl, + C2hhl, + c,h^ + ..) (52) 



where the ellipses denote possible higher order terms |^. In order to investigate 
the relevant degrees of freedom of this theory, we employ again the Stiickelberg 
decomposition, Eq. f|T3|) . We begin by investigating the A9/2 - limit. 



m->0, Mpi -> cx), keeping A9/2 = (m^V'Mp^) fixed. (53) 

Let us for definiteness concentrate on the particular value p = 2 for the coefficient 
in front of the Einstein-Hilbert term, since it corresponds to the case of NMG as 
shown above. Then, in order to recover the FP Lagrangian at the linearized level, 
a = 1/2 should hold. Expanding the Lagrangian of Eq. (!52|) up to cubic order and 
performing the usual rescaling h^^, — )■ h^, + 2ri^i,TT, we obtain, 

A/2 = -l^^''£'^vKl3 



+ -txUtx + {aid^d,nd^'dand''d''7r + a2n7rd^d,TTd''d''7r + (nvr)') , 
^9/2 

where the corresponding coefficients are given as follows, 

tti = -2-2ci-4/3i, (54) 
a2 = 2 - 2c2 + 4/3i - 4/32 + 8/33, 
ag = -2c3 + 4/32 + 8/34. 

In order to avoid the propagation of ghosts as discussed above, we choose these 
coefficients such that the three terms form a total derivative (or vanish) and therefore 
have no effect on the equations of motion. The relevant conditions therefore are, 

2 

«! = "~3'^2 = 203. (55) 



^All terms except the first one in the second line of (15^ can in fact be removed at the 
cubic order by a nonlinear redefinition of the metric perturbation h^^ ; in particular, setting 
h^^ — >■ h^^ + PihP^hpij + [i2hhpy + (i^ih^^haprj^^ + (i^h'^rj^^, will eliminate the hhG^^^-type terms, 
while generating additional /3-dependent polynomial contributions. This leads to an alternative, 
equivalent formulation of the cubic theory; we will however choose to work with the original form 
of the action (l52l) below. 



16 



NMG represents one example of a theory for which these conditions hold. 
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Appendix A 

In this appendix we briefly outline the procedure of integrating out the metric pertur- 
bation hfj^iy in f|T9|) - fl20|) in order to obtain the effective equations of motion, quadratic 
in f^i, (or equivalently, the cubic effective Lagrangian for /^v)- 

For this, we can use the expression for ^G^^ found from f|T9|) in the equation 
for the auxiliary field f l2U]) . and take into account the first-order relation h^y = f^y 
to rewrite everything in terms of f^y only. Specifically, the latter relation is used 
for rewriting the (linearized) Einstein and Ricci tensors in terms of f^y in the first 
line of (fT9|) . as well as the (linearized) quantities such as y^, g^^ and the covariant 
derivatives in the second line of ( IT9|) and in the second term of Eq. f l20lM After 
a rather involved but straightforward computation, one arrives at the following 
effective equation for f^y 

+m' (^Iv'-'r^u - Iv'^f - \rf: + \rf) = o, (a-i) 

where 7^1^ = rf^y + f^y and (a/tC'^'Xt) )*■"'' denotes the contribution to the corre- 
sponding quantity of order n in f^y More explicitly, the first two terms of (]A-ip 
involve the following expressions, 

(v^G"^'^(7))^'^ = ^;;r^ 

^Note that we do not have to go beyond the first order in expressing /i^^ in terms of /^^ for 
obtaining the second-order equations for the latter field. 

^^The perturbative equation of motion for a massless GR graviton h^y is of the form 
{y/gG^"" [9))''^^ + iy/gG'"' {9)Y'^^ + ■■■ = 0. It is crucial that the indices are raised in the latter 
equation; since h^y is defined as the perturbation of the metric g^y (and not its inverse) , the vari- 
ation of the Einstein-Hilbert action w.r.t h^y yields the Einstein's equations with upper indices, 
expressed through the metric perturbation. 
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-□/; - d'^dj) - ^r{d,d^fi: + d^^d^u - - 9^9,/), 

where all indices are contracted with the flat metric. 

It is now straightforward to obtain the cubic Lagrangian which yields (lA-ip as 
the equation of motion for f^^, 



^f2 /■_2f3 _ 1 

fiuJ J fiu 2 



+ i^fU-Vlu-T^f)- (A-II) 



Here, the notation /^^ = f^^f^^ and /^^ = ^'"'fapf^ has been used. 
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